Amplitude death phenomena in delay— coupled Hamiltonian systems 
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Hamiltonian systems, when coupled via time-delayed interactions, do not remain conservative. In 
the uncoupled system, the motion can typically be periodic, quasiperiodic or chaotic. This changes 
drastically when delay coupling is introduced since now attractors can be created in the phase 
space. In particular for sufficiently strong coupling there can be amplitude death (AD), namely 
the stabilization of point attractors and the cessation of oscillatory motion. The approach to the 
state of AD or oscillation death is also accompanied by a phase-flip in the transient dynamics. A 
discussion and analysis of the phenomenology is made through an application to the specific cases 
of harmonic as well as anharmonic coupled oscillators, in particular the Henon-Heiles system. 
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I. INTRODUCTION 

The coupling between dynamical systems can give rise 
to a number of collective phenomena such as synchroniza- 
tion, phase locking, phase shiftin g amplitude death 
[2||, phase-flip 0-Q, hysteresis [l(f, riddling [n| and 
so on [12j. Since communication between the individ- 
ual systems is mediated by signals that can have a finite 
transmission time, many studies account for this by the 
introduction of time-delay in the coupling |13l-[l8|. A 
number of recent studies have examined the manner in 
which delay coupling can affect the collective dynamics, 
particularly since time-delay makes the systems effec- 
tively infinite dimensional [19j . 

When conservative systems are coupled via time- 
delayed interactions, then there are additional consid- 
erations. To start with, the system becomes explicitly 
non-conservative and thus the nature of the dynamics 
changes drastically: in the uncoupled system, the phase 
flow preserves volumes [20j], but in the coupled system 
there can be attractors. This issue is of added interest 
when studying Hamiltonian systems where there can be 
a hierarchy of conserved quantities [2l[ . Studies of cou- 
pled Hamiltonian systems have largely examined the case 
of instantaneous coupling [22j which does not affect the 
Hamiltonian structure. 

In the present work we study time-delay coupled 
Hamiltonian systems and examine the effect of interac- 
tion on the nature of the dynamics. We consider the 
following examples, that of diffusively coupled harmonic 
oscillators that models delay-coupled pendulums, for in- 
stance, and coupled Henon-Heiles oscillators. In the 
absence of coupling, in the former case the motion is 
periodic, while in the latter case the dynamics can be 
(quasi)pcriodic or chaotic. In both instances we find that 
the effect of introducing dissipation is to cause the oscilla- 
tory dynamics to be damped to a fixed point, namely we 
find that there is the so-called amplitude death (AD) Q 
as has been seen in delay-coupled dissipative dynamical 
systems [olfirjl]. 

Although the major effect of the coupling is to make 



the overall system dissipative, there are differences from 
the case when non-Hamiltonian systems are coupled. 
When the dynamics is decaying to a point attractor, there 
is an abrupt transition in the relative phases of the oscil- 
latory transient motion. This is the phase-flip transition 
that has been seen in a number of other systems Q . Here, 
however, there are special values of the time-delay when 
the coupling term effectively vanishes: the underlying 
Hamiltonian structure then becomes apparent. 

Our main results are presented in Sections II and III 
where we consider the cases of coupled harmonic oscilla- 
tors and coupled Hcnon Hcilcs systems respectively. We 
show how AD is reached and the nature of the phase-flip 
transition in both cases. Since the uncoupled systems are 
Hamiltonian, it is possible to define an energy, and while 
this quantity has been studied in coupled feedback oscil- 
lators [23| as a tool to determine onset of AD, its variation 
in the AD regime itself has not been examined. In this 
work we do energy analysis in the AD region and find that 
the energy dissipation is non monotonic as a function of 
the coupling, decaying faster prior to the phase-flip tran- 
sition and slower subsequently. The paper concludes in 
Section IV with a summary and discussion of the results. 



II. DELAY COUPLED HARMONIC 
OSCILLATORS 

The simplest system we consider is that of diffusively 
coupled harmonic oscillators. We consider the following 
equations of motion 



ij + u 2 Xj — e[x k (t — t) — ij(t)] = 



(1) 



where j, fc=l, 2 and j ^ k and Xj and Xj represent the po- 
sition and the velocity of the jth oscillator, and uij is the 
intrinsic frequency. We take the oscillators to be identical 
[24| . u>i = u>2 = co = 1. The parameters e and t repre- 
sent coupling strength and time delay respectively. In 
the absence of delay, e causes the systems to synchronize 
completely. Due to the simple dynamics of the system 
no other significant behavior is observed. When, delay is 
finite the coupling quenches oscillations leading to AD. 



Stability analysis of Eq. ([T]) around the fixed point, 
namely the origin, gives the characteristic equations 



A(A + e)+oj 2 = ±eA exp (-At). 



(2) 



Taking the roots of the Jacobian to be A = a + i/3, the 
condition for marginal stability is a = 0, and substituting 
this condition in Eq. @ we get 

717T nT 

t = t c = — = — - and 

lo 2 

P = ±LJ 

where t c is the critical value at which a = and T is 
the time-period of the uncoupled oscillators. Shown in 
Fig- da) are the first two Lyapunov exponents [25[ of the 
system as a function of r for fixed coupling strength e = 
1. The largest LE is zero only at r c = nT/2 (marked 
in Fig. 1(a) as B and D) and remains negative for all 
other values of r, implying that the system is driven to 
AD except when the delay is an integral multiple of half 
the time period. Further, at these critical delay values 
the system oscillates at the frequency of the uncoupled 
system namely /3 = uj, and the parameter space is di- 
vided into multiple AD regions by the critical delay val- 
ues which are independent of the coupling strength e. 
In contrast, in non-Hamiltonian systems AD islands are 
separated by finite range of delay values @, [l5| where the 
coupling function need not vanish. Hence, in those sys- 
tems the reappearance of oscillations after AD depends 
both on the coupling strength and on the delay, whereas 
in coupled Hamiltonian systems we find that this hap- 
pens only due to delay. 

In the AD regime(s) points of discontinuity in the slope 
of the largest LE (marked in Fig. 1(a) as A and C) indi- 
cate the change in the relative phases of oscillation. This 
is the so-called phase-flip transition Q, and the differ- 
ence in the phases of the coupled subsystems changes by 
7r; see (FigQTj). As in other cases where this phenomenon 
has been observed, there is simultaneously a discontinu- 
ous change in the oscillation frequency [26[, as shown in 

Figntc). 

The uncoupled conservative Hamiltonian systems are 
made dissipativc through the time-delay coupling, and 
this is also reflected in the fact that the sum of all the 
LEs remains negative for all r. Defining the energy of 
the individual oscillators as 



Ej = (xj +njxj)/2 



(3) 



where Q,j is the instantaneous frequency of oscillation, 
the approach to the fixed point can be seen to be at an 
exponential rate in the AD regime 27J as can be seen in 
Fig. [Ha). The exponential decay is however modulated, 
the oscillatory behavior being due to coupling (see the 
inset in the figure). In order to capture the dynamics, we 
define a decay constant as 
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FIG. 1. (Color online) (a) Lyapunov exponents of coupled 
harmonic oscillators Eq. {TJ as a function of the time-delay 
r for fixed e = 1; (b) phase difference between the oscillators 
around point A indicated in (a). Inset figures show trajecto- 
ries xi and X2 as a function of time for the in-phase and the 
out-of-phase dynamics at r =1 and 2 respectively, namely on 
either side of the phase-flip transition, and (c) the frequency 
of oscillation of the subsystems. 



where EJ 1 represents the m-th maxima in the energy time 
series of the j'th oscillator. The averages (-) m and {-)ic 
are performed on m and 100 initial conditions respec- 
tively. 

Since the oscillatory behavior is modulated by expo- 
nential decay, E" l+1 — EJ 1 is also an exponentially decay- 
ing function of m, at rate £j. This rate can be measured 
experimentally and its variation with r is shown in Fig. 
IHh); the variation mirrors that of the frequency change 
at phase-flip, suggesting that energy decays more rapidly 
before the transition than after it. 

At the critical delay values r = nn/uj (see points B 
and D in Fig. HJa)) the largest Lyapunov exponent is 
zero and thus the motion is periodic. Fig. G^a) shows 
orbits for five different initial conditions at the critical 
point; these resemble invariant curves as in conservative 
systems. However, in the vicinity of the critical delay 
values the largest lyapunov exponent is near zero and the 
motion appears periodic after an initial transient. The 
time series of one such periodic orbit is shown in Fig. 
[3fb). Also, as can be seen in the inset, there is very slow 
decrease in amplitude. This occurs since delay is not 
strictly r c . 

Away from delay r c , the dissipation is more pro- 
nounced. The rate of decrease of the amplitude can be 
quantified through the measure 
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FIG. 2. (Color online) (a)Energy of one of the oscillators, 
Eq. ([3|, as a function of time for r = 1. The marked box is 
expanded in the inset, (b) The decay rate, £i as a function of 
r (see Eq. (J2])). Errorbars (marked in red) are calculated for 
100 initial conditions. 
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where xj 1 is the m-th maxima of xj. This is plotted in 
Fig. [3jc) as a function of r in the vicinity of t = 7r/w, 
namely the point B in Fig. [lja). At r — mr/u) the rate 
of decrease of amplitudes approaches zero and hence the 
orbits are almost periodic. Similar behavior is observed 
at point D of FigQJa). 

The reason for reappearance of oscillatory motion is 
straightforward. When the delay is a multiple of half 
the natural period of oscillation, then the coupling term 
effectively vanishes since 

± k (t-T)f»Xj(t) (6) 

and the system effectively becomes conservative. Clearly 
when this occurs, each initial condition gives rise to an 
invariant curve (or, in this case, a nearly invariant curve). 
The better the equality above is realized, the more long 
lived the transients are. 

There is the phase-flip transition at A, C and also 
higher values of r. At each transition a phase differ- 
ence of 7r is introduced between the oscillators resulting 
in anti-phase synchronization at r = T/2 and in-phase 
synchrony as r = T and so on. Hence, the coupling 
function becomes zero at these critical delays. The con- 
secutive oscillatory states alternate between having the 
oscillations being in-phase or out-of-phase, as shown in 
Fig. EJd) (out-of-phase at r = tt/cj, namely at B) and 



FIG. 3. (Color online) (a) Periodic orbits with five different 
initial conditions at r = n/u; (b) time series of position vari- 
able xi as a function of time at r = 3.14. The inset figure 
shows decrease in amplitude of xi; (c) The averaged distance 
Si with errorbars, as function of delay r near r = n/uj. The 
errorbars are calculated for 100 initial conditions; (d) out-of- 
phase motion and (e) in-phase motion in relative phase plane 
x\ — x'2 at r = tt/lu and 2tt/lj respectively (corresponding to 
points B and D of Fig. [H a ))- 

Fig. [3^e) (in-phase at r = 2ir/ui at the point marked 
D) respectively. Note that the phase-relation between 
oscillators is independent of the initial conditions. 



III. COUPLED HENON-HEILES SYSTEMS 

In order to examine the dynamics when the uncoupled 
systems are capable of exhibiting chaotic motion, we ex- 
amine the behavior of two non-integrable Henon-Hciles 
systems [2fj| . 

ij + Xj + IxjDj - e(±k{r — t) — ±j[t)) = 

jjj + Vj + x) - y) = (7) 

As is well-known, in the uncoupled case (e=0) the 
system has both regular and irregular behavior largely 
depending on the total energy as well as on the initial 
condition [2(| . Shown in Fig. 2] are the Poincare maps 
for two different initial conditions, one leading to regular 




FIG. 4. (Color online) (a) Poincare map of two orbits in 
the Henon-Heiles system, Eq. @, at energy Ej = 0.13. Red 
(outer) and black (inner) dots represent regular and chaotic 
motion respectively. 

motion (red points), while another leads to chaotic dy- 
namics (black points), at the same energy Ej = 0.13 just 
below the dissociation limit Ej = 1/6. 

A. Instantaneous coupling (r = 0) 

In the absence of delay r = 0, Eq. (JT)) reduces to 
the case of simple diffusive coupling. The effect of in- 
creasing the coupling strength, namely e, is to induce 
simplicity to the resulting collective dynamics. Shown in 
Fig.© are the the fraction of initial conditions (/) lead- 
ing to quasiperiodic motion. We take 100 pairs of ran- 
dom initial conditions from the bounded region of phase 
space (Fig. @|. In one case, when the initial condition 
pairs of quasiperiodic motions are taken, the collective 
dynamics due to interaction remains quasiperiodic (solid 
line (black circles): QP + QP). However, if the initial 
motion is chaotic then the resulting dynamics becomes 
quasiperiodic only after certain value of coupling strength 
(dashed line (red triangles): C + C). Similar behavior 
is observed when initial conditions of mixed chaotic and 
quasiperiodic motions are used (dotted line (green stars): 
C + QP). These results indicate that for small coupling 
strength the Hamiltonian structure still exists, but for 
larger values of coupling the collective dynamics becomes 
quasiperiodic; in this sense coupling induces simplicity in 
such coupled systems. 



B. Time delay coupling 

When the two systems arc coupled in presence of delay, 
r (initial conditions taken from cither the regular or irreg- 
ular motion) the largest Lyapunov exponent quickly be- 
comes negative. This is shown as a function of the delay 
in Fig. [BJa): almost as soon as the delay is switched on 
the oscillators are driven to AD state. The phenomenol- 
ogy of this higher dimensional system is very similar to 
that of the coupled harmonic oscillators: the largest LE 
has the same characteristic shape — there is a change in 
the slope at point A, where it is clear that the phase- 
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FIG. 5. (Color online) Fraction of initial conditions that lead 
to quasiperiodic motion. The solid (black circles), dashed 
(red triangles) and dotted (green stars) curves are the results 
for the different cases of initial conditions corresponding to 
quasiperiodic (QP + QP), chaotic (C + C) and mixed chaotic 
and quasiperiodic (C + QP) motions respectively. Averages 
have been taken over a sample of 100 initial conditions. 



flip transition occurs. Shown in Fig. [SJb) is the com- 
mon frequency of oscillations as a function of delay which 
changes discontinuously at r ~1.65. Since the phase of 
the oscillators is not clearly defined in such systems, we 
infer the phase relation from the time-series (inset fig- 
ures) of the two systems in the neighborhood of the point 
of discontinuity. The phase difference changes from to 
7T along with the frequency jump as in the simpler 1— 
dimensional harmonic system. 

Here also energy decreases exponentially in this AD 
region. We define the energy of individual systems in the 
usual manner (20| 

_ xl + yl xl + yl , y) 
&3 ~ 2 2 ' ' —■ W 

and quantify the energy dissipation as in Eq. © by a 
quantifier £j . The variation of this quantity with r can 
be seen in Fig. [6jc). It confirms that energy dissipation 
is faster before phase flip transition and slower thereafter. 

When the largest Lyapunov exponents approach zero 
(at points marked B and D in Fig. GHa)) the motion is 
oscillatory, decaying very slowly to the fixed point. Crit- 
ical delays are at half the natural time period of oscilla- 
tion, and since the natural frequencies of the oscillators 
is equal to T', the time period is 2ir. The Poincare sec- 
tions of representative trajectories at t — 3.14 arc shown 
in Fig. [TJa) which arc quasiperiodic (cf. Section HA). In 
the vicinity of this point, the rate of decay can be com- 
puted numerically as discussed earlier, and the quantifier 
8, Eq. ([5]) defined to locate the oscillatory state also has 
a minimum for at r m 3.14. At successive critical points 
(t = 3.14 and t = 6.28) the quasiperiodic motions alter- 
nate in the nature of the phase synchrony; see Fig.[7Jd,c). 
Unlike the case of harmonic oscillators, though, the cou- 
pling term does not quite vanish, so the emergence of 
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FIG. 6. (Color online) (a) Largest Lyapunov exponent as 
a function of the time-delay, r, and (b) the oscillation fre- 
quency, for coupling strength e = 1. The inset shows trajec- 
tories xi(t) and X2(t) for r =1.4 and 1.8 respectively, (c) The 
decay constant £i as a function of the delay r. 



oscillations is not as pronounced in this case. 



IV. SUMMARY AND DISCUSSION 

In the present paper we have explored the effects of 
incorporating time-delay in the coupling between Hamil- 
tonian systems. The coupled system effectively becomes 
dissipative, and in the absence of other attractors, the 
total system displays amplitude death. However the cou- 
pling can effectively vanish for specific values of the time- 
delay, at which delay the the system naturally appears to 
be conservative. We find that this behaviour differs from 
that of delay-coupled non-Hamiltonian systems, where 
AD islands are separated by finite ranges of delay values 
where the coupling function need not vanish. 



Orbits near these critical delay values reflect both dissi- 
pative and conservative behavior: different initial condi- 
tions give rise to different (seemingly) invariant curves 
which are very long lived transients. Energy dissipa- 
tion in the AD regime is found to be associated with 
the phase flip transition, and the damping is faster prior 
to the transition, and slower after it. 

The finite velocity with which signals are transmitted 
gives rise to intrinsic delays in the coupling, and as such 
this is germane in both dissipative as well as conserva- 
tive systems. Nevertheless, the effect of delay has been 
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FIG. 7. (Color online) (a) Poincare section at r = 3.14 
for three different initial conditions; (b) time series of xi; (c) 
variation of Si as a function of delay. The error bars are 
computed from a large number of initial conditions, (d) The 
relative phase in xi — X2 plane to show anti-phase motion at 
t — n and (e) in-phase motion at 2tt. 

explored to a limited extent in conservative systems. The 
results presented here may have more general applicabil- 
ity in coupled systems with other conservation laws such 
as in ecological contexts [1^] . 
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